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We consider here the percolation problem in thin films, both in the direction normal to the film 
and in the direction parallel to the film. We thereby describe here the cross-over between 2D and 3D 
percolation, which we do on cubic and square lattices. The main relations are derived using scaling 
and real space renormalisation arguments. They are checked by numerical simulations, which also 
provide the numerical prefactors. We calculate in particular the correlation length parallel to the 
film, the average mass and the mass distribution n(m) of the clusters. In particular, we show that 
the latter is given by a master function of h~ D+1 ^ a2l ' a \p — PcCi)| " 2 m, where h is the thickness of 
the film and D, 1/3, 02 are tabulated 2D and 3D critical exponents. p c (h) is the percolation threshold 
of the film which we also calculate. These results are of interest in particular for describing the glass 
transition in thin polymer films. 
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I. INTRODUCTION 



The building of macroscopic continuous objects (ag- 
gregates, or clusters) by the random dispersion of parti- 
cles (sites) or bonds has been the subject of many studies 
over the past decades, and has been formalised by various 
models of percolation 0, H H B The corresponding is- 
sues are essential for determining important macroscopic 
properties such as the electric conductivity Q, the 
visco-elastic behaviour : 5| of various systems, or more 
generally mechanical properties |(J . For a review on the 
latter aspect, see e.g. Q. The main quantities of inter- 
est have been expressed as power laws, close to a critical 
point, and the main results are summarised by a list of 
critical exponents, which do not depend on the details of 
the models but only on the spatial dimension. Due to the 
dependence of critical phenomena on spatial dimension, 
the cross-over between 2D and 3D behaviours has been 
the subject of various studies, e.g. in the case of mag- 
netic phenomena for describing the phase transition 
behavior of thin magnetic films for instance. For studying 
the conductivity of thin composite systems, Clerc et al 
considered the cross-over between 2D and 3D percolation 
and calculated the percolation threshold as a function of 
the thickness of the film, using finite scaling arguments 
0- This issue has also been considered by Vicsek and 
Kertesz [T3 - 

Another field where percolation appears to be a key 
concept is the glass transition in supercooled liquids. 
The most prominent feature of this phenomena is a dra- 
matic increase of the viscosity when cooling such liquids. 
Though very steep, this increase is continuous and does 
not appear to involve a phase transition. For a review 
on the glass transition see e.g. 0,^2- To account for 
experiments which have demonstrated the heterogeneous 
nature of the dynamics when approachin g th e glass tran- 
sition [3|ll EE 13 E3 ( see e-g- El Ell for review), 



it has been proposed recently 21, 22] that the glass tran- 
sition corresponds to the percolation of domains with re- 
laxation times larger or equal to the arbitrary time r g set 
for defining the glass transition. According to this model, 
the slow domains correspond to denser regions resulting 
from density fluctuations of Gaussian statistics. These 
slow domains, of typical size 2 to 4 nm, coexist with 
faster domains with relaxation times order of magnitude 
shorter 

El EillllUli^. This model > aimed at explain- 
ing experimental results regarding the heterogeneous na- 
ture of the dynamics close to the glass transition, allows 
also to explain the shift of the glass transition tempera- 
ture T g of thin polymer films, a few tens of nanometer 
thick. Indeed, it has been demonstrated over the past 
ten years that the dynamical behaviour of thin polymer 
films is very different from that of the same polymer in 
the bulk. It is now well established that films with weak 
interactions with their susbstrate (or freely suspended 
films) display a reduction of T„, of about 20 K for films 



10 nm thick |23l l24j. On the contrary films with strong 
interactions with their substrate display an increase of 
T g by as much as 60-K" for films 10 nm thick |25l]2q.|27J. 
For review of the corresponding issues see e.g. |28l l29j| . 
In the model of Long and co-workers [2l], 123. l30| , these 
effects result from percolation mechanisms. According to 
this model, the glass transition is controlled by the perco- 
lation of small subunits of relaxation time r g . Therefore, 
the shifts in T g in a film are closely related to the dif- 
ferences in the percolation properties in a film (i.e. in a 
system of small, finite thickness) with respect to a bulk 
system. In the case of a thin suspended film, percolating 
in the direction parallel to the film requires a larger frac- 
tion of these slow subunits, and thereby takes place at a 
temperature lower than the bulk T g . 

In the case of a strongly interacting film, the glass tran- 
sition occurs when aggregates of slow subunits have a 
diameter comparable to the thickness of the film, due 
to different boundary conditions. As a consequence, in 
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both cases the glass transition corresponds to the tem- 
perature at which the correlation length of the 3D perco- 
lation problem is equal (or comparable) to the thickness 
of the film, in one case (weak interactions) above the 3D 
percolation threshold, and in the other one (strong in- 
teractions) below the 3D percolation threshold. In the 
case of intermediate interactions between the monomers 
and the substrate 0, |13 , as described in [3(j , the glass 
transition corresponds to a lateral extension of a slow ag- 
gregate larger than the film thickness, and such that the 
number of monomers of this aggregate in contact with 
the substrate is large enough so that the time for this 
aggregate to desorb is equal to r g . Describing all these 
cases -weak interactions, strong interactions, intermedi- 
ate interactions- requires thereby a precise understan ding 
of the cross-over between 2D and 3D percolation |2lll22| . 
This is the aim of this paper. A better description of this 
cross-over should also be useful for describing the me- 
chanical properties of films made of composite materials, 
or for the conductivity properties of thin films Q, . 

The questions which we address here stem from our 
modeling of the glass transition in thin films. When con- 
sidering a thin suspended film, we argued 0, that 
the glass transition corresponds to the percolation in the 
direction parallel to the film. Then, we will consider this 
problem here, in films of dimcnsionlcss thickness h, ex- 
pressed in units of the size £ ~ 2-4 nm of the dynamical 
heterogeneities 0, l^- Then, one of the questions is: 
what is the percolation threshold p c (h) knowing the 3D 
percolation threshold p^. D l How does the lateral exten- 
sion of the aggregates increase when approaching this 
threshold? When considering films strongly attached to 
their substrate, we argued that the glass transition cor- 
responds to the appearence of a fraction of order 1 of 
monomers from one interface connected to the opposite 
interface by continuous path of slow subunits. We thus 
study the correlation function between sites on one inter- 
face and sites on the other, to determine precisely how 
the connected fraction evolves with the fraction p of oc- 
cupied sites. Finally, to consider the crossover between 
these two regimes, we need to know the size of the ag- 
gregates, the number of monomers of one aggregate in 
contact with the interacting substrate, and the distribu- 
tion of mass of the aggregates. 

The present work is based on the idea that a system 
of finite thickness h may be renormalised to a 2D sys- 
tem, by changing the 2D percolation threshold p 2 c D into 
a renormalised value p c (h). A first step will be therefore 
to determine p c (h). We describe a real-space renormal- 
isation procedure to transform the quasi-2D problem of 
a film of finite thickness into a 2D problem. This allows 
for calculating various quantities of interest, as a function 
of known 2D and 3D critical exponents and of a renor- 
malised occupation probability. The corresponding re- 
sults are exact at the level of scaling laws. To check these 
predictions, and to obtain the values of the pre-factors, 
we then perform numerical simulations. The paper is or- 
ganized as follows. In section II. A, we introduce some 



TABLE I: The values of the percolation thresholds in various 
cases 



Lattice 


Site 


Bond 


square (2D) 


0.592746 


0.50 


simple cubic 


0.3117 


0.2492 


FCC 


0.198 


0.119 



classical definitions, notations and we recall some basic 
concepts and results regarding percolation theory, cither 
in 2D or in 3D. Specifically we recall the various quan- 
tities which have a critical behaviour at the percolation 
threshold. In section II. B, we introduce some classical 
finite size scaling arguments regarding the definition of 
the percolation threshold on a finite system. Then in sec- 
tion II. C, we show how the percolation problem in a film 
of finite thickness can be mapped on a 2D system with 
a new percolation threshold p c {h). Various scaling laws 
are derived, regarding p c (h), the mass of the aggregates, 
the connectivity between both interfaces. In section III 
we describe the numerical Monte-Carlo algorithm used 
for our simulations. Finally in sections IV, V and VI 
we discuss the results of our numerical simulations, in 
connection with the various scaling laws derived in II. C. 



II. BACKGROUND ON PERCOLATION 
A. Critical behaviour at p — p c 

Let us recall some basics about 2D or 3D percolation 
which will be useful in the following. For a review of the 
definitions and results presented here see e.g. Q. The 
control parameter in site percolation problems is the site 
occupation probability p. The percolation threshold is 
the critical value p c at which an infinite cluster appears. 
The value p c is not universal. It depends on the partic- 
ular type of lattice which is considered. Values of the 
percolation thresholds in various cases are summarized 
in Tabic [I] In the following, the various exponents are 
generically designed by greek letters. When some rela- 
tions are valid only using a 2D or a 3D exponent, the 
corresponding exponents are written with a subscript 2 
or 3 according to the dimension. 

For the sake of definiteness and when it is necessary to 
be specific, we will consider in this paper percolation on 
cubic lattices (3D or in films), and on square lattices for 
2D systems. Only the prefactors depend on this choice, 
the exponents of the various scaling laws that we derive 
do not. We call n(m) the number of clusters of mass m 
(the mass m of a cluster is the number of sites which 
belong to it) per lattice site (i.e. the total number of 
clusters of size m divided by the volume of the system) , at 
a given occupation probability p. The quantity mn{m) 
therefore equals the probability for a given site to belong 
to a cluster of mass m. The occupation probability is 
thus p = Y]mn(m). The distribution of cluster masses 
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is represented by a function of the form: 
n(m) ~ 7ti t f(m/m^) 
where f(m/m^) behaves like 

/(m/m c ) ss e- m ' m < 



(1) 



(2) 



for to 3> to^ and has a finite value at zero. For to 
smaller than the typical value ttiq, the power law domi- 
nates. For m larger than m^, the exponential dominates, 
which means that clusters with a mass larger than rriQ 
are exponentially rare, diverges at p c with a critical 
exponent denoted by 



™C w \p - p c 



-l/a 



(3) 



In 2D, the exponent <r 2 = 36/91 = 0.3956 (or l/a 2 = 
2.5278) and the exponent r 2 = 187/91 = 2.055 (see Ta- 
ble P). 

On the other hand, a mean cluster mass may be defined 
for p < p c as the second moment: 



M = 



m 2 n(m) 
y~] m n(m) 



(4) 



The mean cluster mass M diverges at the threshold p 
p c with another exponent 7: 



M ^\p-pc 



(5) 



In 2D: 72 = 43/18 = 2.389. 

Since ^ran(m) = p (or equivalently, since r > 2), the 
divergence is contained in the numerator in eq. (@J , which 
reads: 



M 



J m 2 - T f (jj^j dm 



Note that M grows with \p — p c | more slowly than . 
This is because M is an average mass (second moment 
of the cluster number n(m)), whereas is the cut-off 
mass of the distribution. Combining eq. 10, JSJ and 
one finds: 
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(7) 



Note that, at p = p c , the distribution n(m) reduces to 
n(m) ss mT T . 

We consider now the spatial extension of the clusters. 
One can define the correlation function (or connectiv- 
ity function) g(r) as the probability that a site at the 
distance r from an occupied site is also occupied and be- 
longs to the same cluster. The correlation length may be 
defined from g(r) as: 



(8) 



For p < pc, g(r) is represented by a function of the form: 



9(r) 



-{d-2+rj) 



G(r/0 



(9) 



where d denotes the dimensionality of space, and 77 is the 
anomalous exponent of the correlation function [l| . For 
r < (, the power law dominates. This means that at 
short scale, the system is insensitive to the presence of a 
cut-off length at a larger scale. At large scale (r ^> £), the 
function G decays exponentially G(r/() ss exp(— r/(). 
The correlation length £, which is proportional to some 
typical cluster diameter, diverges at the threshold as: 



C ~ \P ~ Pc 



(10) 



where, in 2D: 1/2 = 4/3, and in 3D: ^3 rj 0.88. 

The radius of a cluster of mass m may be defined as: 



2R„ 



E- 



(11) 



where 



is the squared distance between two occupied 



sites which belong to the same cluster of mass m, and 
where the sum is then averaged other all clusters of mass 
to. According to the definitions, it is easy to check that 
one has: 



c 2 



^2 to 2 n(m) 



(12) 



The correlation function g(r) is related to the structure 
of the clusters, i.e. to the distribution of the mass inside 
a cluster. In 3-D, let us draw a volume of size h (of 
volume ft 3 ) within a large cluster, such that h is small 
with respect to £. The number of sites within this volume 
which belong to this cluster is: 



(6) And for h < (: 



g(r)r 2 dr 



(13) 



(14) 



In 3D, the exponent rj is relatively small: 77 = —0.068. 
Setting h = £ in eq. I|14|) . one obtains the average mass 
M, as defined in eq. @: M« ( 2 ~ v , which gives another 
relation between critical exponents: 



(2 -r])v = 7 



(15) 



We consider now the distribution of radius of the clus- 
ters for p < p c . The radius R m of a cluster is related to 
its mass to by a scaling law: 



R„ 



(16) 



However, the exponent (or apparent exponent) D is not 
the same below, at and above p c 0- For p — p c , the 3D 
critical exponent is D3 « 2.53, and for p < p c , the scaling 
relation, i.e. the fractal exponent, depends on the con- 
sidered radius. Thus, more generally, it is assumed that 
the relation between the mass and the spatial extension 
is of the form 0] : 



Rn 



m^ D f[(p-pc)m° 



(17) 
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TABLE II: The values of some of the critical exponents, which 
will be used in the following 



Exponent d — 2 d — 3 

~ 187/91=2.055 2l8 

a 36/91=0.3956 0.45 

7 43/18=2.389 1.80 

v 4/3 0.88 

l/D(p = p c ) 49/91=0.5275 1/2.53=0.3952 

D(p<p c ) 1.56 2 

rj 0.2084 -0.068 



where the function / has to be determined numerically 
and satisfies the following properties: as m < rj 
\p — p c \~ 1 / (T (where ttlq has been defined in eq. 10), i.e. 
asi= \p— p c \m a -C 1, /(a;) approaches a constant value. 
For x » 1, /(x) ~ x 1 /- '" 1 /- , where D' is the fractal 
exponent of large clusters for p < p c . Away from p c , 
m C ~ \p — Pc| is finite, and defines a crossover from 
the behaviour R m rj w}I d for to < rrif to another power 
law i? m rj to 1 / 15 for to > TO^. At p c , to^ diverges. The 
exponent l/£) in cq. 117(1 is therefore the one which is 
measured at p c (or very close to p c ). This exponent takes 
the values l/D 3 = 1/2.53 = 0.39526 in 3D and 1/D 2 = 
48/91 = 0.5275 in 2D. For clusters with a mass larger 
than TO£, at p < p c , the value of the fractal exponent 
is smaller: D$ = 2. This means that for p < p c , large 
clusters are more tenuous than small ones, while they are 
self-similar whatever the size at p c . The above discussion 
means that the crossover takes place on a spatial scale 
(. For small clusters with R m < £ (or for small masses 
to < mf), the value of the critical exponent is that at 
p = p c (D rj 2.53 in 3D), and for larger clusters the 
exponent is that at p < p c (D = 2 in 3D). The cross-over 
between these two behaviours takes place at R m = £j 
corresponding to the mass to^ = C, D , where D is the 
critical exponent at p c . Note that the cross-over goes 
to larger scales when p tends to p c and there is a single 
exponent at p c . All this means that at short scales (i.e. 
at scales significantly shorter than the correlation length 
C), the system does not know whether it percolates at 
larger scales or not. 

Combining eqs. J3J), (|10|l and (|16fl . one gets another 
relation between critical exponents: 

- = vD (18) 

(T 

where D has the value at p = p c here. 



B. Percolation in a finite system 

The problem in numerical simulations is: how to ex- 
tract the value of the percolation threshold (which is de- 
fined at infinite size) from simulations done in boxes of 
finite size ? Percolation in a system of size L is defined by 



the appearance of one cluster of radius L. The probabil- 
ity that a system of linear size L (volume L 3 ) percolates 
at an occupation probability p is of the form: 

*l(p)« f(L(p c ~p) v ) (19) 

where the exponent v has been introduced in eq. Ill) I I. 
Eq. (|19fl means that the probability to percolate in a 
system of size L is identical to the probability of appear- 
ance of a cluster of radius L in an infinite system, or, 
equivalently, that a system of size L percolates when the 
correlation length = L (where £ is defined e.g. in 
cq. JTT3J))- According to eq. JT5J, the curves &l(p) be- 
come steeper as L increases. We call p c {L) the occupa- 
tion probability corresponding to a given probability to 
percolate P. It follows from eq. (|19|l that in a system of 
size L, Pc(L) varies as: 

\ Pc {L)-p c \ rjL- 1 /" (20) 

Eq. (|20|l will be used in the following to determine per- 
colation thresholds p c (extrapolated at infinite size) from 
simulations in boxes of various finite sizes. The percola- 
tion threshold p c (n) will be given by the number giving 
the best regression with exponent —1/v (according to 
eq. (J23). 

C. Rescaling a system of finite thickness h 

Consider a thin film with mesh size unity, probability 
of site occupation p and of thickness h. The probability 
of occupation p may be larger or smaller than p 3D and is 
such that the corresponding 3D correlation length is 
larger than the thickness h of the film, that is: 

Pa(h)<p< P t(h) (21) 

The quantities p£(h) and p^(h) correspond to a correla- 
tion length for the 3D percolation problem comparable 
to the thickness of the film, is smaller than the 3D 
percolation threshold p 3 c D , and p+ is larger than p 3 c D . We 
aim now at calculating the percolation threshold of the 
film, p c {h). Note that we consider here percolation in 
the direction parallel to the film. By a coarse graining 
procedure, this problem can be transformed into a 2D 
problem. The corresponding coarse graining function is 
denoted by $. After one iteration, the probability of oc- 
cupation is transformed into <I>(p). p^P is the fixed point 
of <I>. Let us define K > 1 as the scale factor of the renor- 
malisation procedure. TV is a number such that 

k n = h (22) 

which means that the film is transformed into a 2D 
squared lattice after N iterations. If the initial film is at 
the percolation threshold, one requires that after these 
N iterations the film is transformed into a 2D lattice 
close to the 2D percolation threshold, that is p = p c (h) 
is transformed into & N (p c (h)) such that 

& N \p c {h))-pl D ~pl D -pl D (23) 
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FIG. 1: Schematics illustrating the renormalization of the 
percolation problem in a film of thickness h to a two dimen- 
sional system. 



One has then 



\ N ip c (h)-pi u ) = a 



where a is a positive number comparable to p 2D 
and where 

X-—( n -n 3D ) 



(24) 
(25) 



Therefore 



Pc(h) « P l u + a\ 



-N 



P, 



:\D 



+ ah~ x ' V3 (26) 



where z/3 is the 3D critical exponent for the correlation 
length, which satisfies the relation: 



V3 



logK 

log A 



(27) 



Note that the film percolation threshold p c (h) is larger 
than p 3D -percolating in thin films is more difficult than 
in the bulk- and that at p c (h) the 3D correlation length is 
comparable to h. In the following we take p£(h) = p c {h). 
Let us consider now a situation where the probability of 
site occupation in the film is p < p c (h). By the same 
argument, p is transformed into <$>( N> *(p) with 



„3D 



\ n (p-pV j ) 



3D\ 



(28) 



This procedure is illustrated in Fig. ^ One has then 

& N \ Pc (h))-& N \p) W X N ( Pc (h)-p) « h^ipcity-p) 

(29) 

The quasi-2D problem of the film with probability of oc- 
cupation p on a lattice of mesh size unity is then mapped 
on a 2D lattice with probability of occupation $W(p) 
and mesh size h. The average mass M'{h) of the ag- 
gregates, expressed in numbers of supersites of size h, is 



therefore: 

M'(h) « Mj|#W(p c (/0) - * W b)r 72 (30) 
wM^^lp-pcWl)" 7 " (31) 

where M' is a prefactor of order unity. Note that the 
above result insures that for p — p~(h), i.e. when the 
3D correlation length is equal to the thickness of the 
film, the typical mass of the aggregates on the scale h is 
of order one as it should be. Note indeed that p~ (h) is 
given by a relation 



Pc ( h ) = P 



3D 



(32) 



where a' is a positive number of order 1. Let us study 
now the structure of the clusters in the film, below but 
close to the percolation threshold p c (h), more precisely 
in the situation h <C C3-D- In this case, the lateral exten- 
sion of the dominating cluster is large with respect to h. 
The 2D correlation length is, expressed in the equivalent 
renormalised 2D system: 

h}' v *\p-p c {h)\) (33) 

Coming back to the original system, in units of elemen- 
tary sites: 



Cn«/i(ft 1/l,3 b-PcWl)" 



(34) 



Here £|| is the correlation length in the direction parallel 
to the film and provides the dominant lateral extension 
of the clusters. Of course this result is meaningful only 
when the occupation probability p is such that the 3D 
correlation length is larger than the thickness of the film. 
In units of the elementary sites of the initial film, M' 
should then be replaced by the average mass within a 
/i-cube. This is done using eq. ifTijl: Mq = M Q h 2 ~ r > (be- 
cause in all the interval studied here h < C3d)- Mq is 
a number of order unity. Finally, the average mass per 
cluster is, in units of elementary sites: 



M 



M Q h 2 ^(h}/ V3 \p-p c {h)^ 72 (35) 



Let us consider now the distribution of clusters n(m). 
It obeys a scaling relation similar to eq. (|17|l (see also 
eq. (QJ). In the equivalent renormalised 2D system, the 
number of clusters of mass m' (in units of /i-cubes) is (see 
eqs. Q and (SJl): 



n(m') w m T2 g2 



(36) 



with: 



a- 2 

(37) 
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cj2 is the 2D function as defined by equation (f), and has 
to be determined numerically. r 2 = 187/91 = 2.055 and 
l/o- 2 = 91/36 = 2.5278 are the 2D values of the expo- 
nents (see Table [J. The expressions for the masses M' 
and M can be obtained in a different way as that used 
above, which we consider now. From the mass distribu- 
tion n(m'), one deduces the relation 



M' a mf~ T2 



(38) 



Using the same argument as that used to obtain equation 
(|3~5)) . one obtains 



M m h 2 - m mf~ T2 



(39) 



Then, from equation l|37|) , the average mass M' on scale 
h can be written as 



M 1 a h-^-^/^^lp - Pc (h)\- (3 - T2 ^ c 



(40) 



Finally using the relation (3 — r<i)j(Ji = 72, one recovers 



M' a ^ 72/l/3 |p-p c (/i)r 72 



(41) 



which is equation (|3 1|) . To come back to the original film 
of thickness h, we need now to obtain the distribution 
n(m) of the mass m of the aggregates on the scale of 
elementary sites. When the system is rescaled by a factor 
k (in 3D), the correlation length is rescaled as C i- 
Since the relation = must be preserved, 

is thus rescaled as: 

TO£ 



^1/(0-3^3) 



(42) 



With l/{av) = D (sec cq. (TL%1) V this leads to the rescal- 
ing rriQ = h D3 m' ( -, where D 3 is here the 3D value of the 
exponent at p = p c . Thus, we have to consider a scaling 
relation of the type: 



(to) » h UJ m- T2 g 2D [/i- jD3+1 /( ff2l/3 )|p-p c (^)| 1 / CT2 TO 



or equivalcntly 



n(m) ^ h^m' T2 g 2D (h- D ^) 
where u> can be determined by writing: 



M 



m 2 n(m)dm a h^+^^m' 3 ^ 2 



(43) 
(44) 

(45) 



Comparing with eq. (|39|l . one obtains the equation u> + 
Ds(3 — T2) — 2 — 1J3. Thereby the distribution of cluster 
masses is expected to be described by equation (|4*3"|) or 
equation (|44l) , the value of the exponent u> being given 
by 



2-V3-D 3 (3-r 2 ) 



(46) 



Thus, the quantity hr^m T2 n(jn) should give a master 
curve when plotted versus the quantity h^ D3+1 ^' 72 ^ \p— 
p c (/i)| 1 / <T2 m. The exponent u> has the value lu = —0.323. 
In the following, we consider numerical simulations re- 
garding the relations derived in this section. 



III. DESCRIPTION OF THE SYSTEM AND 
THE ALGORITHM 



We compute the site percolation threshold in a square 
(2D) or simple cubic (3D) lattice of variable linear dimen- 
sions. Unit sites are occupied at random with a probabil- 
ity p. An improved random number generator (uniform 
deviate provided in numerical recipes in Fortran) , is used, 
which, thanks to a shuffling procedure of the output, is 
known to provide perfect random numbers, within the 
limits of the floating point precision. Indeed, the ran- 
dom generator has to be called a number of times equal 
to the volume of the system times the number of sample 
systems used to average the data (that is, at most a few 
10 11 ), which precludes the use of a standard, low quality 
random number generator. 

The algorithm used to determine percolation thresh- 
olds is itself based on a renormalization procedure, which 
is either two- or three-dimensional. For instance, at a 
given step in the 3D procedure, d-cubes of linear size 
d = 2 m (of volume V — d 3 ) are generated. The faces 
of the cube are occupied by atomes belonging to vari- 
ous clusters. Then, 8 such (i-cubes are gathered to form 
one cube of linear size 2d = 2 m+1 (volume 8V), in which 
the clusters from all 8 elementary cubes are connected to- 
gether and re-numbered. A Hoshen-Kopelman algorithm 
is used to renumber the clusters [33| . Then only the new 
clusters connected to one external face of the new cube 
at least are considered for the next step. The initial step 
consists in connecting together 8 unit sites of the lattice. 
Iterating the 3D procedure to times gives a system (a 
cube) of size 2 m . Then, a film of thickness h and lateral 
dimension d = h p may be generated by iterating p times 
a 2D procedure similar to the 3D one. The system per- 
colates when there exists at least one cluster connecting 
the 6 external faces. The probability to percolate P{p) is 
averaged over a number of different configurations (typ- 
ically 50 to 500). Using the renormalization procedure, 
the curves P(p) for all values of d (of the form d = 2 J ) 
are obtained in a single run. Simulations were run on 
Pentium IV Personal Computers operating at 1.4 GHz 
with 512 Mo RAM or on a Dec Alpha Work Station. 

The algorithm was validated by studying the well- 
known cases of 2D (resp. 3D) square (resp. simple cubic) 
systems, in order to verify that percolation thresholds are 
obtained with a satisfactory precision. The probability 
to percolate P is measured as a function of the occupa- 
tion probability p. In two dimensions, the curves P(p) 
obtained for various sizes of the system (from d = 2 3 = 8 
to d — 2 16 = 65536, which corresponds to a surface d 2 
from 64 to 4.295 x 10 9 ) become steeper when d increases, 
in agreement with eq. H19|) . We find here that all curves 
intersect at a fixed point which corresponds to an oc- 
cupation probability p 2D « 0.5927 and a probability to 
percolate P{p) ~ 0.327. The value p c is in agreement 
with that given in 0. The curves P(p) are well fitted 
(at least in the vicinity of the percolation threshold) by 
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a function of the form: 



P(JP) = \ 



1 + tanh 



P-Pc(d) 
9(d) 



(47) 



in which both adjustable parameters p c (d) and q(d) de- 
pend on the size d of the system. The form of this func- 
tion is purely heuristic, as it merely provides a unam- 
biguous way to extract the value of the threshold for 
each value of d. When the quantity \og(p c (d) — pj, D ) 
is plotted as a function of logd, a straight line with a 
slope —0.75 ± 0.01 is obtained. This is in concordance 
with eq. (J2DJ, with the theoretical value —l/v 2 = —3/4. 
The quantity q(d), plotted as a function of d in log- 
arithmic scale does also give the power law behaviour 
q(d) w d^ 1 ^ 2 , in agreement with eq. fTT)|) . In 3D as 
well as in 2D, the curves P(j>) obtained for various sizes 
of the system (from d = 8 to d = 2 9 = 512, i.e. a 
volume 512 to 1.342 x 10 8 ) strongly depend on d. All 
curves intersect at the fixed point p z c D — 0.3117 and a 
probability to percolate P(p) ~ 0.073. The value pl° 
is in agreement with that given in [lj. The curves P{p) 
are fitted by a function of the same form as in the 2D 
case (eq. 1^3) • In 3D as well, the quantity p c (d) — pl° 
varies as d^ 1 ^ 3 , with the expected 3D value of the ex- 
ponent I/1/3 = 1/0.88 = 1.1364. Note that the prefactor 
in this curve (as well as in the 2D curve) is arbitrary, as 
it depends on the probability level P at which p c (d) is 
measured. The scaling behaviour q{d) ~ d^ 1 ^ 3 is well 
verified in this case also, in accordance with eq. I|19|) . 



IV. PERCOLATION THRESHOLD IN THE 
DIRECTION PARALLEL TO A FILM OF 
THICKNESS h 

We consider the variation of the percolation threshold 
in a system of dimensions d x d x h as a function of the 
thickness h of the system. Note that the corresponding 
results are of interest for calcula ting the negative shift 
of T g in a freely suspended films [2l|, |3(| • The thickness 
h was varied from h = 1 (which corresponds to the 2D 
system) to h — 128. For each value of h, the dimension 
d was varied from d = h to d = 1024. The algorithm 
combines 3D renormalization up to a cube of dimensions 
h x h x h, then 2D renormalization of /i-cubes up to a 
system of dimensions dxdxh. Percolation is defined here 
when one cluster at least connects all 6 external faces of 
the volume dxdxh (see Fig. 0) . For each couple of values 
(d,h), the probability to percolate P is measured as a 
function of the occupation probability p. Some examples 
of the curves P(p) are shown in Fig. [3] For each value 
of h, the curves corresponding to different values of d 
intersect at a fixed point, which may be considered as 
the percolation threshold p c (h). The curves are fitted by 
the same heuristic function as before, with two adjustable 
parameters p c {d, h) and q(d, h) : 



Pip) = \ 



1 + tanh 



p - p c (d, h) 
q(d, h) 



d 



h 



m 




FIG. 2: Percolation in a system of dimensions dxdxh: there 
is a cluster spanning the whole surface of the system. This 
figure illustrates the case of a non-adsorbing film. 
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(48) 



0.315 0.320 0.325 

occupation probability p 

FIG. 3: The probability to percolate P as a function of the 
occupation probability p, obtained in systems of thickness h 
(h — 32 and h — 64) and lateral extension d: □ : d = 1024, 
o : d = 512, o : d = 256, A : d = 128 ,V : d = 64, + : d = 32. 
For each value of h, the curves P(p) become steeper as d 
increases and intersect at a fixed point p c (h). 



The number value giving the best linear regression 
p c {d,h) - p c (h) w d- 1 /" 2 (h being fixed), with the 2D 
value of the exponent v 2 = 4/3, will be adopted for the 
percolation threshold p c (h) in a film of thickness h, in 
the limit d — > oo. Within error bars, this value p c {h) 
coincides with the fixed point mentionned above in the 
ensemble of curves P(p). Thus, for a given thickness h 
and varying lateral size d, the following variation is ob- 
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TABLE III: The values of the percolation thresholds in a film 
of thickness ft 



h p c {h) 

2 0.47424 

4 0.3997 

8 0.3557 

16 0.3329 

32 0.3219 

64 0.3165 




log h 

& io 

FIG. 4: ■ : the prefactor fx(h) in the curves p c {h, d) — p c (h) 
versus d, as a function of ft, o : the prefactor go (ft) in the ad- 
justable parameter q(d,h), in logarithmic scale. The straight 
lines show the predicted behaviour fj,(h),qo(h) oc ft -0,386 . 

served: 

Pc (d,h)~ Pc (h)^ ^d- 1 ^ 2 (49) 

with i>2 = 4/3. The values obtained for the percolation 
thresholds p c {h) in a film of thickness ft are summarized 
in Table IIIII The variation of the prefactor /j,(h) may 
be estimated as follows. The 3D behaviour (up to the 
size ft) gives p c (d) -p\ D = Ad- 1 /" 3 (see eq. ®). Thus, 
setting d — h in eq. (|49|) : 

Pc {h) + ^{h)h- 1 ' 1 ' 2 « p 3 c D + Ah- 1 / 1 " 3 (50) 

which gives, since p^P — p c {h) « ft, -1 / 1 ' 3 : 

fx(h) » h l ' V2 - l ' Va (51) 

The curve p,(h) is shown in Fig. 0] ■ The variation is 
compatible with the exponent 1/^2 — 1/^3 = —0.386. 
Note that the prefactor in eq. (|5l"1) is arbitrary, since it 
depends on the probability level P which is adopted to 
define percolation. 

The same argument may be used to study the varia- 
tion of the fitting parameter q(d, ft). It is indeed observed 




0.0 0.5 1.0 1.5 2.0 

log h 

& io 

FIG. 5: The reduced percolation threshold p c (ft) — p^P as 
a function of the thickness ft of the system, in logarithmic 
scale. The straight line represents the scaling law with the 
3D exponent —I/V3 = —1.1364. The scaling law is observed 
asymptotically. 



that, varying d for a fixed value of ft, a 2D scaling be- 
haviour is observed: q(d,h) = qo{h)d~ 1 / v ' 2 (see Fig. 0J. 
The prefactor is observed to vary as qo(h) = qoh~ Q 3S6 , 
in accordance with the expected variation, which may be 
infered by exactly the same argument as was used for 
/i(ft) (eq. H51|) h The prefactor qo is not arbitrary here: it 
has the value 0.565. In Fig. O the percolation threshold 
Pc(h) — Pc D i s plotted as a function of ft in logarithmic 
scale. The best linear fit has a slope —1.1364 which corre- 
sponds to the 3D value 1/^3 = 1/0.88 of the exponent, in 
agreement with eq. (12611 . Thus, in the asymptotic limit, 
the curve p c (h) writes : 

Pc (h)-pl D *ah- 1138 (52) 

From Fig. [21 the prefactor a has the value a « 0.45. Note 
that this value is of the same order as, but different from 
the value p 2 c D - pl° = 0.281 (cf eq. Q). 

V. CORRELATION FUNCTION BETWEEN 
BOTH INTERFACES OF A FILM OF 
THICKNESS h 

We consider here the correlation function between both 
interfaces of a thin film. The regime of interest here is 
that of a thickness h larger than the 3D correlation length 
of the percolation problem. Note that the corresponding 
results are of interest for calcula ting the positive shift of 
T g in a strongly adsorbed film [2l], y(J]. Indeed, in the 
case of a strongly interacting film, the glass transition 
is related to the presence of one percolating cluster con- 
necting the upper surface to the lower one (see Fig.[SJ). In 
the sense of percolation, glass transition will thus occur 
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1.0- 



0.5- 



FIG. 6: Schematics of a film of thickness h in the regime £ < 
h. Only percolation clusters connected to the lower surface 
(i.e. to the substrate) are shown. Very few clusters reach the 
upper surface. 



o 



0.0- 



in this case, as soon as the 3D correlation length £ be- 
comes of the order of (still being smaller than) of h. The 
effective percolation threshold defined in this way is thus 
shifted downwards with respect to p^. D (since p^. D corre- 
sponds to £ — > oo). Let us consider the average number 
of sites to which an occupied site of a given interface is 
connected to the other interface. This number n c (h) can 
be calculated by integrating the (3D) correlation func- 
tion, expressed in eq. ©: 



n c (h) 



i+na ('- 2 +ft 2 ) 



rdr 

p-" 3 e-<dp (53) 



Integrating by parts, this gives: 



n c (h) 



<dp (54) 



The second term is negligible for h 3> C- 

In our simulations we have calculated the number 
n s (h) of sites of one interface which are connected to the 
other interface with at least one site. We expect a similar 
dependence of n c (h) as that given for n s (h) by equation 
eq. (154(1 . In our simulations, we consider a system of di- 
mension d x d x h (with d large with respect to h). For 
a given value of the occupation probability p such that 
£ < h, we compute the number n s (h). The value of d is 
d = 1024, h ranges from h = 1 to h = 64. The algorithm 
makes use of a ID renormalisation procedure : 2D slices 
of size d x d are generated and then glued together re- 
cursively, leading to systems of thickness 2, then 2 2 , etc, 
up to 2 6 = 64. 

For a given value of the occupation probability p (with 
p < Pc°), the quantity log n a {h) exhibits a linear be- 
haviour in the limit of large h values. Each curve was fit- 
ted with the first term in eq. (|54|) , with two independent 
adjustable parameters 773 and £. The values obtained for 
C in this way differ by less than 15% from those obtained 
by simply taking the slope of the asymptotic linear vari- 
ation at large h values, i.e. by neglecting the influence of 
the h~ ,n factor (remember that the exponent 77 is quite 
small). This difference is systematic and does not alter 
the variation as a function of p^. D — p as it is shown in 
Fig. 7. 



-0.5- 



-2.0 



-1.5 
log p 

& 10 J c 



-1.0 



3D. 



FIG. 7: The correlation length £(p) obtained from a linear 
fit of the curves n s (h) at large h values (see eq. 1541 '). ver- 



sus Pc — p, in logarithmic scale (p^ 
threshold p 3 c D =0.3117). 



is the 3D percolation 



In Fig. [7| the correlation length £(p) is plotted as a 
function of p^. D — p, in logarithmic scale. The data are 
compared to a straight line of slope —vj, = —0.88, which 
is the behaviour expected in the limit p — > p c . An excel- 
lent agreement is obtained in the asymptotic limit. Thus 
we have: 



c(p)«Co(pr-p) 



(55) 



with Co « 0.185. Note that two opposite requirements 
have to be reconciled in this section: observing the crit- 
ical power law £(p) « (p c — p)~" 3 urges for coming close 
enough to p^P , while maintaining £ small enough to make 
the regime h > £ observable. 



VI. CLUSTER NUMBER IN A FILM OF 
THICKNESS h 

The regime of interest here corresponds to a 3D corre- 
lation length larger than the thickness of the film, either 
below the 3D percolation threshold, or above. The dia- 
gram in Fig. [S] summarizes the results obtained so far. p 
is the occupation probability, h the film thickness (that 
is, the number of slices in the simulated system). The 
curve Pc{h) is defined as the ensemble of points where 
the measured 3D correlation length £ is equal to h. Note 
that this curve corresponds to the percolation threshold 
associated to glass transition in a strongly adsorbed film, 
as it was defined in Section According to eq. 1551 it is 
described by the equation: 



„3£> 



h 



-1/1/3 



(56) 
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FIG. 8: Diagram of the percolation thresholds Pc(h) and 
pt{h) as a function of the film thickness h. Pc{h) is the per- 
colation threshold in an adsorbed film. It is defined here as 
the ensemble of points for which C 3 °(p) = h ■ Diamonds 
symbols correspond to the same data set as in Fig. [7] The 
continuous line is the fit with equation 1561 pj{h) is the per- 
colation threshold in a freely suspended film. Square symbols 
corresponds to the data in Fig [5] The dashed line is the fit 
with eq. 1521 in section IIVI 

The curve p£(h) is denned by eq. It corresponds to 
the percolation threshold in the direction parallel to the 
film. Note that it corresponds also to the glass transition 
in a freely suspended film of thickness h, as defined in 
Section llVl The two curves Pc(h) and pt{h) define 3 re- 
gions in the plane (p, h). In region I, one has C < h < d. 
Clusters connecting both faces of a film of thickness h 
are exponentially rare. In Region III, there is an infinite 
cluster extending throughout a film of thickness h. Re- 
gion II corresponds to intermediate cases h < £ < d. For 
a given thickness h, this corresponds to going from p~£ (h) 
to Pc(h) within region II. In this region, some clusters 
have a lateral extension large with respect to h. As one 
comes closer to the curve p+(h), the correlation length 
parallel to the plane of the film, i.e. the lateral extension 
of clusters £ diverges. One can see such an aggregate in 
Fig. |5J Thus, in this section, we study the distribution 
of cluster sizes and masses, both on the surface and in 
the bulk of a film of thickness h. The corresponding re- 
sults are of interest for describing the cross-over between 
freely suspended films and strongly adsorbed films as far 
as glass transition is concerned |30j. 

We will check that the distribution of clusters and their 
structure may be described by the procedure detailed in 
Section fll CI in which we assumed that a percolation sys- 
tem of thickness h may be renormalized to a true 2D 
system by coarse-graining at the scale h. Note that the 
system has to be larger than the typical lateral extension 
of large clusters, which become larger as one comes closer 
to p+(h). On the other hand, the second moment of the 



h 




































































































































































































































































































































































































































































































































































































































































































FIG. 9: Schematics of a cluster percolating through the film 
thickness, in the regime h < £ < d. The 'surface mass' s of 
the cluster is given by the ensemble of sites at the surface of 
the film, which are emphasized. 

distribution of cluster masses n(s) is dominated by the 
contribution from large clusters, which are exponentially 
rare. Thus, we need here to simulate both large enough 
systems and a large enough number of realisations of a 
given system, which is quite demanding in terms of com- 
puter size and time. The volume V = h x d 2 of the 
simulated systems is of the order a few 10 s typically. Re- 
sults are averaged over 20 to 50 different realisations, 
depending on the distance to the p^(h) curve. 

A. The average mass of the clusters 

Let us consider a cluster in a film of thickness h, in the 
regime h < £ < d (i.e. in Region II of the diagram in 
Fig.|HJ. It was shown in Section 111 CI how a film of thick- 
ness h may be renormalized to a 2D percolation system. 
Specifically, in the regime h < £ < d, the average mass 
M, expressed in units of elementary sites and defined as 
in eq. Q, is given by eq. (|33)l . i.e.: 

M(h,p) = M h 2 -" 3 (V^b-PeWl)" 72 ( 57) 

The 2D behaviour is reflected in the exponent —72. We 
have obtained the distribution of cluster masses n(m) for 
different values of p and h (within region II) in systems 
of total volume of the order 2.6 x 10 s , and computed the 
second moment M(h,p) (see eq. (QJ). 

Let us first consider the variation versus p, h being 
fixed. For each value of h, the variation of M(h,p) 
versus p c {h) — p gives a power law with the exponent 
—72 = —2.389 in the limit p close to p c (h), in agreement 
with the expected behaviour. To illustrate this power 
law and the influence of the thickness h of the system, 
the curves M(h,p) are plotted as a function of the quan- 
tity X = h^fjPcQi) -p) in Fig. HHi for different val- 
ues of h. The values at \ogX = for each portion of 
curve in Fig. 1101 (corresponding to a given value of h) 
gives the prefactor Mo(h) in the curves M(p,h) versus 
X. In the limit of large h, the behaviour of the prefactor 
Mo{h) is compatible with a power law Mq(K) = Moh 2 ~ m 
(2 — 773 = 2.068). The prefactor Mq has the value 
M w 0.233. 

The quantity M' = h- 2+V3 M(h,p) plotted versus X 
should then give a master curve of the form M' — 
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FIG. 10: a: the average mass M(h,p) vs X — 
h 1/u3 (p c (h) -p), b: the quantity M' = h- 2+V3 M(h,p) vs 
X, for different values of h: □ : h — 1, o : h — 2, o : h — 4, 
® : h = 8, A : h — 16, V : h = 32. The straight lines corre- 
sponds to the exponent —72 = 43/18 = 2.389. In b, a master 
curve is obtained to an excellent approximation. 



M' Q X~^. This is illustrated in Fig. ITUb. The expected 
behaviour is indeed observed to an excellent approxima- 
tion in a large range of variation (more than 4 decades in 
M) . The full scaling behaviour of M(h, p) as a function of 
p c {h) —p and h (eq. (57)) is thus summarized in Fig. llOb . 
Remember that M' represents the mass expressed in 
number of renormalizcd supersites. Alternatively, the 
same work may be done with the surface mass s of the 
clusters, i.e. the number of sites belonging to a cluster 
which are in contact with the substrate (the surface mass 
of clusters is denoted by s in order to distinguish it from 
their total mass m. See Fig. EJ. The average surface mass 
was obtained in the following way. The distribution n(s), 
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FIG. 11: The reduced average surface mass h 1+r,s S(h,p) vs 
the quantity X — h 1 '" 3 (p c (h) — p) in log scale, for different 
values of the thickness h. The straight line corresponds to the 
exponent -72 = 43/18 = 2.389. 



which gives the number of clusters with a surface mass 
s, was recorded for different values of h and p (in region 
II of the diagram in Fig. The average surface mass 
S(p, h) was then computed as the second moment: 



S(h,p) 



EjMf) 
E sn ( s ) 



(58) 



where the summation extends over all clusters which per- 
colate through the film thickness. 

In a way similar to M(h,p), S{h,p) is expected to scale 
like: 



S{h,p) = S h 1 -^ {h l ' v *\p-p c {h)\) 



-72 



(59) 



When plotted versus (p — p c (h)), for each value of h, 
S(h,p) gives a power law with the exponent —72 = 
—2.389, in agreement with the expected behaviour and 
with the behaviour observed for the quantity M(p,h). 
From eq. (|59|l . a master curve should be obtained 
when the quantity h 713 ' 1 S(h,p) is plotted versus X = 
fi 1 /" 3 (p c (h) —p). This is shown in Fig. ^] Here also, 
the expected asymptotic behaviour is observed to an ex- 
cellent approximation in a large range of variation for S. 
The prefactor 5o has the value So rs 0.24. 



B. The distribution of cluster masses in the surface 
of the film 

In this Section, we show that not only the average 
mass, but also the whole distribution of cluster masses 
obeys scaling laws in the two variables h and p — p c {h), 
which allows to draw a master curve. According to 
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Eq. H43|) . the number of clusters of mass to obeys (in 
the limit of large m): 



eq. ( I60|l , obeys the following scaling law: 



h u m T2 n{m) w g 2 D 



,-D 3 + l/a 2 v 3 



\p-Pc(h)\ 1/a *m 



(60) 

i.e., the quantity h "m T2 n(m) should give a master 
curve when plotted versus the quantity h~ D3+1 ' a21/3 \p — 
p c (/i)| 1 / <T2 m. The characteristic exponents of h have the 
values -D 3 + \ja 2 v z = 0.342 and ui = -0.323. Con- 
sider first the variation versus p — p c (h), h being fixed. 
The distributions of clusters, represented by the quan- 
tity m T2 n(m), obtained for h = 4 and different values of 
p —pc(h), are plotted as a function of to in Fig. 112b. For 
each value of h [h = 1 to h = 32), an excellent superpo- 
sition property is observed when m T2 n(m) is plotted as a 
function of the reduced variable \p — p c (h)\ 1 /' J2 m. This is 
illustrated in Fig. 112b . in which the various curves shown 
in Fig. IT2*k (h — 4) are superposed. 

Two exponents have to be adjusted in order to su- 
perpose the ensembles of curves obtained for different 
values of h. First, plotting the quantity log m T2 n{m) 
as a function of the reduced variable fi^ L>3+1 ' c ' 2U3 \p — 
p c (/i)| 1 / <T2 m allows the slopes obtained in the asymp- 
totic limit at large to to coincide. Then, the quantity 
log h~ u 'm T2 n(m) is plotted as a function of the reduced 
variable h~ D3+x l a2V3 \p - p c [h)\^^ 2 m. This procedure is 
illustrated in Fig. 112b . As mentionned above, according 
to eq. HJn| a master curve should be obtained. 

The quantity log h^^rn T2 n{m) does indeed show a lin- 
ear asymptotic behaviour at large m values, i.e. for clus- 
ters larger than the correlation length Q. It is observed 
in Fig. 112*1 that the slopes obtained at large m values 
indeed coincide, within numerical uncertainties. In con- 
trast, vertical superposition of the curves obtained for 
various h values is only observed asymptotically, i.e. to- 
wards the largest investigated value h — 32. The large 
dispersion in the curves comes from the numerical lim- 
itations of the present simulations, especially for large 
values of h and/or p c (h) — p. In this regime, where the 
correlation length ( is large, the lateral size d of the sys- 
tem has to be large, and the distributions n(m) extend 
towards large m values. As already mentionned, com- 
puting the distributions n(m) with a good sampling of 
large to values is very demanding in terms of memory 
size and computer time. In any case, it can be concluded 
at this point that our data are compatible with the scal- 
ing law in eq. 1601 in the asymptotic regime of large h 
values, which is particularly clear when considering the 
data corresponding to the thicknesses h = 8, 16 and 32 
respectively. 

To illustrate further the scaling properties expressed 
in eq.( l6T"|) . we have fitted the exponential decrease of 
m T2 n(m) at large m values. This procedure gives the 
quantity m ( ^(p,h) 1 which depends on p and h and may 
be considered also as a typical mass describing the distri- 
bution of the clusters within the film (see eq.([""J)). From 



-1/0-2 



(61) 



First, for each value of h, it is observed that m^(p,h) 
follows a scaling law when plotted as a function of 
{p — Pc{h)), with the exponent — l/oa = —2.528. Then 
a master curve is obtained when h~ D3 rriQ(p,h) is plot- 
ted as a function of the quantity h 1 ^ 3 (p — p c (h)). The 
obtained master curve is shown in Fig. 1131 



C. The structure of the clusters in a film of 
thickness h 

The structure of the clusters is described by the rela- 
tionship between the mass m of a cluster (the number 
of sites belonging to the cluster) and its size r. Clusters 
have a fractal structure characterized by the exponent 
D defined as r w m x / D , in which the value of D effec- 
tively measured varies according to the location of the 
system with respect to p c , or equivalently, according to 
the typical range of sizes which is probed compared to 
the correlation length ((p) (see eq. ifTTjl and the discus- 
sion thereafter) . The cluster size may be defined as either 
the gyration radius or some measure of the cluster exten- 
sion. 

We have studied the relationship between the size and 
mass of clusters in a film of thickness h, within the inter- 
val Pc(h) < P < pt{h). The squared size of a cluster i is 
defined here as r 2 = Axj + Ay 2 , where Axi and Aj/i are 
the maximum extensions of the cluster i along x and y 
(x and y are the coordinates parallel to the plane of the 
film) , or more precisely, the maximum extensions of the 
intersection of the cluster with the limiting surfaces of 
the film. The mass rrii is the total number of sites which 
belong to cluster i (within the whole volume of the film) . 
The quantity rf is then averaged over all clusters of mass 
m, in order to obtain the squared average size r 2 (m) as 
a function of m: 



r 2 (m) 



1 



Vn(m) 



E 5 



(62) 



where the sum is extended over the Vn(m) clusters hav- 
ing the mass m (V = d 2 x h is the volume of the system). 
The relationship between the mass to and the average 
size R = (r 2 (m)) 1 / 2 of the clusters in a film of variable 
thickness h is illustrated in Fig. 1141 

Some examples of the curves i?(m), obtained in a sys- 
tem of lateral dimension d = 2048 and different values 
of the thickness h (from h = 1 to h = 16) are plotted in 
logarithmic scale in Fig. 1141 for an occupation probabil- 
ity p — 0.315. The observations which can be made on 
Fig-El are the following. First, two different regimes are 
observed. They are characterised by well defined power 
laws, with different exponents. At small R values (typi- 
cally, R smaller than or of the order of h), the mass grows 
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FIG. 12: Scaling properties of the distribution m T2 n(m), 
where n(m) is the number of clusters of mass m per lat- 
tice site . a: log m T n(m) plotted versus m for a thickness 
ft = 4 and different values of the occupation probability p: 
□ : p = 0.34, o : p = 0.36, A : p = 0.365. b: log m T n(m) 
plotted versus \p — p c (ft)| 1//<T2 m, for different values of ft : A 
: ft = 1, o and ♦ : ft = 4, □ : ft = 16. In this represen- 
tation, the curves obtained for different values of p at fixed 
ft superpose, c: ft - " log m T n(m) plotted versus the quan- 
\p — p c (h)\ 1 ^' T2 m, for different values of p and ft: 
1, o : p = 0.34, ft = 4, o : p = 0.34, ft = 8, 
= 16. V : n = 0.31. ft = 32. The curves 



tity ft' 
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FIG. 13: The quantity obtained by fitting the distribution 
m T2 n(m) at large m. The quantity ft _D3 m^ is plotted versus 
X = h 1/v3 (p-p c (ft)), for different values of ft: □ : ft = 1, 
o : ft = 2, o : ft = 4, A : ft = 8, V : ft = 16. A master curve 
is obtained to an excellent approximation. The straight line 
corresDonds to the exDonent — 1/ao = —2.528. 
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FIG. 14: The average lateral size R = [Ax 2 + AY 2 

plotted as a function of the total mass m of the clusters, in 
films of various thickness ft: □ : ft = 1, o : ft = 2, o : ft = 4, 
A : ft = 8, V : ft = 16. The occupation probability p — 0.315. 
The straight lines have a slope 1/1.561, which corresponds to 
the fractal exponent Di = 1.56 at p < p c . 



with R faster than for larger R values. Indeed, clusters 
should behave in a 3D way in this regime. Within the 
whole interval pj(ft) < p < pt(h)> t nc 3D correlation 
length is larger than ft. Thus, the expected behaviour in 
this regime should be described by the fractal exponent 
I/-D3 at Pc D , according to eq. I|17|l and the discussion 
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thereafter. This is indeed the case: the slope of this part 
of the curves is compatible with the value I/-D3 = 0.395 
at p = p 3 c D . This regime is illustrated by the dashed line 
in Fig. Then, at R « h, each curve has a crossover 
to another regime. For R > h, one enters a quasi 2D 
regime, in which the correlation length in the equivalent 
(renormalized) 2D system is small. Therefore, the slope 
in this regime has to be compared to the value of the 
fractal exponent at p < p 2 c D : 1/Df = 0.641 Q. This 
quasi 2D regime is illustrated by the solid line in Fig. ITU 
The observed behaviour is indeed compatible with the 
exponent l/D^- 

Thus, the results in Fig. E] illustrates the quasi 2D 
behaviour of the film: a film of thickness h behaves as 
a strictly 2D film. In the equivalent renormalized 2D 
system, the 2D scaling law, eq. I16l writes: 

rri = m R' D * (63) 

where too is a prefactor of order unity (the superscript 
< indicates that the exponent D% is that at p < p 2D ). 
Coming back to the original system of thickness h, the 
size (expressed in units of elementary sites) becomes 
R' — > R = hR' and the mass is transformed as to' — > 
to = h D3 m', since £30 > h. Eq. I|63|) thus leads to the 
following scaling: 

/R\ D i 

to = m h D3 ( — J (64) 

From the data in Fig. O the value of the prefactor 
too w 0.22. Close inspection of the data in Fig. 1141 shows 
however that the apparent exponent in the regime R > h 
is slightly smaller than l/Df in some range of R val- 
ues at the beginning of this regime. Indeed, the scal- 
ing to' w i?' D 2 is only valid when R > 0i , where 0i is 
the correlation length parallel to the film. In the range 
h < R < 0\ , one should observe a power law to cx i? 1//£>2 
where I/D2 = 49/91 w 0.538 is the critical exponent at 
p c . It just happens here that the range is too small for 
this regime to be clearly visible. For instance, 0i (ex- 
pressed in units of /i-cubes) is: 

C|| «Co(h 1 ^\p+(h)-p\y Va (65) 



which gives for the particular case of the data shown in 
Fig- El a value Qi « 3.2 (with a prefactor £ of the order 
1). 



VII. CONCLUSION 



We have described here the cross-over between 2D and 
3D percolation. This work demonstrates that a system 
of finite thickness h may be renormalized to a strictly 
2D system for describing the percolation properties, pro- 
vided that the 2D percolation threshold p 2D is renor- 
malized to an appropriate value p c (h). The values p c {h) 
have been determined numerically with a good precision 
for h = 2 to h = 128 (see Table lllljl . in the case of a 
cubic lattice. We have studied 3 different regimes, which 
describe 3 different physical situations. These regimes 
are represented in Fig. At a given value of the occu- 
pation probability p, if the thickness h is larger than the 
3D correlation length, the problem is a true 3D problem. 
The fraction of sites of one interface connected to the 
other interface by a cluster is exponentially small. The 
corresponding correlation function and correlation length 
are given by 3D critical exponents. Our numerical sim- 
ulations provide the corresponding pre- factors. If, at a 
given value of p, the thickness h is smaller than the 3D 
correlation length, then the problem is equivalent to a 2D 
problem. We provide the corresponding transformations 
for calculating the correlation length parallel to the film, 
the average mass of the clusters and the cluster mass dis- 
tribution, among other relations. These quantities can 
be expressed as functions of known 2D and 3D critical 
exponents. Our numerical simulations provide the corre- 
sponding pre-factors for these quantities, as well as the 
master function for the mass distribution. These results 
can be used for describing t he g lass transition in thin 
polymer films as described in y(j , and should be also of 
interest for describing mechanical or electrical properties 
of films made of composite and disordered materials. 
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